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O ' Abstract 

We investigate the observational signatures of three models of the early Universe in the -B-mode polar- 
ization of the Cosmic Microwave Background (CMB) radiation. In addition to the standard single field 
inflationary model, we also consider the constraints obtainable on the loop quantum cosmology model (from 
Loop Quantum Gravity) and on cosmic strings, expected to be copiously produced during the latter stages 



of Brane inflation. We first examine the observational features of the three models, and then use cur- 



fvq , rent .B-mode polarization data from the BICEP and QUaD experiments to constrain their parameters. 

£> ' 

We also examine the detectability of the primordial B-mode signal predicted by these models and fore- 



CO 

these three models of the early Universe predict different features in the CMB B-mode polarization power 



cast the parameter constraints achievable with future CMB polarization experiments. We find that: (a) 



spectrum, which are potentially distinguishable from the CMB experiments; (b) since B-mode polarization 
measurements arc mostly unaffected by parameter degeneracies, they provide the cleanest probe of these 
early Universe models; (c) using the BICEP and QUaD data we obtain the following parameter constraints: 



r = 0.021q 26 (lcr for the tensor-to-scalar ratio in the single field inflationary model); m < 1.36 x 10~ 8 M p i 
and fc* < 2.43 x 10 _4 Mpc _1 (lcr for the mass and scale parameters in the loop quantum cosmology model); 
and G[i < 5.77 x 10~ 7 (lcr for the cosmic string tension); (d) future CMB observations (both satellite mis- 
sions and forthcoming sub-orbital experiments) will provide much more rigorous tests of these early Universe 
models. 
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I. INTRODUCTION 

Observations of the Cosmic Microwave Background (CMB) radiation have proved a valuable 
tool for studying the physics of the very early Universe. Scalar, vector and tensor perturbations 
generated in the early Universe have left observable imprints in the temperature and polarization 
anisotropies of the CMB. Recent experiments, including the Wilkinson Microwave Anisotropy 



Probe (WMAP) satellite |lHg|, QUaD |5(, BICEP [6j] and others [7J4ll|], have led to a precise 
determination of the basic parameters of the standard ACDM cosmological model, including the 
parameters describing the primordial density perturbations. 

According to this concordance model, the Universe underwent a period of near-exponential 
expansion, termed inflation, at very early times. The standard model of inflation is based on the 
single field slow-roll scenario. In this scenario, the expansion is driven by a scalar field (the inflaton) 
gradually rolling down a flat potential during the inflationary stage. Inflation ended when the slow- 
roll conditions were broken, and the inflaton decayed into relativistic particles which re-heated the 
Universe. 

In spite of many phenomenological successes of inflation based on effective field theory, serious 
problems remain concerning the origin of the scalar field driving inflation, namely the singularity 
problem [12] and the trans-Plankian problem [13]. Consequently, efforts have been made to realize 
inflation in a more natural way from some fundamental theory of microscopic physics. Brane 
inflation [14], [l5[] from high dimensional string theory is a typical example. In this scenario, the 
Universe is embedded into a high dimensional warped space-time. The anti-brane is fixed at the 
bottom of a warped throat, while the brane is mobile and experiences a small attractive force 
towards the anti-brane. Inflation ends when the brane and the anti-brane collide and annihilate, 

ision, cosmic strings would be copiously 



16, 



171 ] . Searching for this cosmic string 



initiating the hot big bang epoch. During the brane colli 

produced, and would leave an imprint on the CMB sky 

signal in the CMB is an important way to test the correctness of this scenario. 

Another approach for realizin g a period of inflation, based on Loop Quantum Gravity (LQG) 



has been proposed recently (see 18| for instance). LQG is a non-perturbative and background- 
independent quantization of General Relativity. Based on a canonical approach, it uses Ashtekar 
variables, namely SU(2) valued connections and conjugate densitized triads. The quantization 
is obtained through holonomies of the connections and fluxes of the densitized triads. More im- 
portantly, when the energy density of the Universe was approaching the critical density p c , the 
Universe entered into a bouncing period due to repulsive quantum geometrical effects. Thus, the 



big bang is replaced by a "big bounce". This, to some extent, avoids the singularity problem in 
the standard ACDM model. 

Differentiating between these three classes of models (standard field inflation, brane inflation 
and loop quantum cosmology), which are motivated by different microscopic physics, is a crucially 
important goal for modern cosmology. Since the primordial scalar, vector and tensor perturbations 
produced in these models are quite different from each other, they will in general leave different 
signatures in the CMB radiation. Numerous authors have constrained the parameters of standard 
9 eld i»«o» m „ d e ls from CMB and large SC a,e _ „,— ,. (S ee for e Xampl e Q, Q, 
20] for some recent analyses.) Most of these analyses have made use of both temperature and 
polarization CMB measurements and their constraints have been dominated by the temperature 
measurements. However, with the advent of a new generation of CMB polarization measurements 
[a, Ifij], it is now possible to obtain meaningful constraints solely from measurements of the B- 
mode polarization of the CMB. Since .B-mode polarization on very large scale is generated only by 
tensor perturbations in the early Universe, this is a particularly attractive technique: a detection 
of S-mode polarization at large enough angular scales must be due to gravitational waves (tensor 
perturbations) and the connection with early Universe physics is then very clear. In addition, 
small scale S-mode polarization is possibly sourced by cosmic string. This potential of .B-mode 
measurements to constrain inflation parameters was demonstrated recently to great effect by the 
BICEP collaboration who used their 5-mode constraints on large angular scales to obtain a 95% 
upper limit on the tensor-to-scalar ratio of r < 0.73. Note that on small scales, l?-modes are also 
generated by gravitational lensing of the dominant .E-mode polarization signal and so in general, 
large-angular scale measurements are required in order to avoid confusion from the signal due to 
lensing. 

In this paper, we extend the investigation of inflationary constraints from CMB B-mode po- 
larization measurements alone by considering the constraints obtainable on the three classes of 
models described above. The paper is organized as follows. In Section [II Al we briefly review the 
characterization of CMB polarization in terms of E- and i3-modes. In Section III Bl we present 
the currently available B-mode constraints, and the predicted noise levels for some current and 
future CMB experiments. In Section IIIQ we first present the likelihood and hyper-parameter 
analysis methods, which we will use in the parameter estimation. We then describe the Fisher 
information matrix formalism which we will use to forecast the constraints obtainable using future 
experiments. In Section IIII Al we discuss the tensor perturbations which arise in the single field 



inflationary (hereafter SFI) model 



831 ] , and in Section IIII F3|, we constrain the parameters of the 



SFI model using the BICEP and QUaD data. In Section IIIICl we calculate the single-to-noise 
ratio for a number of future CMB experiments and their combinations, and present forecasts for 
the constraints obtainable on the tensor-to-scalar ratio r with future observations. In Section nil Dl 
we discuss four types of single field slow-roll inflation models, and their detectability with future 
experiments. We then follow a similar line of discussion for the LQG model in Section HVl and the 
brane inflation/cosmic string model in Section [V] We summarize our results in Section IVIi 

II. 5-MODE POLARIZATION AND ITS OBSERVATIONS 

A. B-mode polarization 

Let us first briefly review the statistics of the CMB polarization field. The polarized radiation 



field can be described by a 2 x 2 intensity matrix iy(n) [211 ]. where n denotes the direction on the 
sky, and Iij(n) is defined with respect to the orthogonal basis (ei, e2) which is perpendicular to n. 
Linear polarization is related to the two Stokes parameters, Q = 4(111-/22) and U = 2-^121 whereas 
the temperature anisotropy is T = j{I\i + ^22)- The polarization magnitude and orientation are 
given by P = ^/Q 2 + U 2 and a = \t&rT l {U/Q). 

As spin ±2 fields, the Stokes parameters Q and U change under a rotation by an angle V as 
(Q±iU)(n) — > e Zf2l 'v(Q±iU)(n). Thus, (Q±iU)(ri) requires an expansion with spin ±2 spherical 



harmonics 



22] 



lm 

The multipole coefficients a) can be calculated as 



(Q±zU)(n) = Y,4t 2) [±2Y lm (n)\. (1) 



«/™ 2) = J(Q ± ^)(n)[ ±2 l^(n)]dn. (2) 

The E- and S-mode multipoles are defined in terms of the coefficients a) in the following manner: 

„^ _ _I r „(2) , „(-2), B _ _]_ ( (2) _ (-2), () 

a im- ^\i mlm '' lm ~ 2i lm lm 
One can now define the electric polarization sky map E(n) and the magnetic polarization sky map 
.B(n) as 

E(n) = Y,af m Y lm (n), S(n) =^a^ m (n). (4) 

lm lm 

The scalar field E(n) and the pseudoscalar field B(n) completely describe the polarization 
field. E'-modes are curl-free modes and appear as symmetric radial and tangential polarization 
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FIG. 1: Comparison of different theoretical predictions, and the currently available data for the B-mode 
power spectrum. The unit for m is M p i and fc* is Mpc" 1 . 

patterns on the sky. S-modes are divergence-free modes with left-handed and right-handed vortical 
polarization patterns on the sky. 

One constructs the various CMB power spectra by correlating the T, E and B modes in harmonic 
space. In the absence of parity- violating effects 8J|, there are only four non-zero cross-correlations: 



TT, TE, EE and BB. The angular power spectra of the polarization fields are defined as 



21 



where the brackets denote an ensemble average. If the fluctuations are Gaussian distributed, all of 
the cosmological information is encoded in the angular power spectra of Eq. (jSJ). 



B. Constraints on the i?-mode signal 

In Fig. [H we plot the current constraints on the B-mode power spectrum along with some rep- 
resentative B-mode signals from different theories. The red curve is the predicted C BB calculated 
from LQG, with a mass parameter m = 10 -8 M p i and k* = 0.002 Mpc" 1 (see Section [IV Al for a 
discussion on LQG). The blue curve is the predicted C BB from SFI for a tensor-to-scalar ratio 
r = 0.03, and the pink curve is the BB power spectrum generated by cosmic strings for a string 
tension Gfi = 10~ 7 and wiggling parameter a = 1.9 (see Section |V| for a discussion on cosmic 
strings). The green curve is the B-mode signal from gravitational lensing which acts as a source 
of confusion when attempting to measure the primordial .B-mode signal. The black, purple and 
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FIG. 2: Polarization noise power spectra for forthcoming experiments. Note that these curves include 
uncertainties associated with the instrumental beam. The blue curves show the -B-mode power spectrum 
for the standard inflationary model with r = 0.03. In the left panel, we plot the instrumental noise for 
the Planck and CMBPol satellites, as well as the lensing .B-mode signal and the noise level for the ideal 
experiment. In the right panel, we plot the instrumental noise for the ground-based PolarBear and QUIET 
experiments, and the balloon-borne Spider experiment. 



orange points with associated error-bars are the currently available -B-mode data from the WMAP 
5-year observations (/ < 20 [23]), the BICEP experiment (9 band powers [a, I24J), and the QUaD 



experiment (23 band powers 5 



25]). 



The WMAP constraints are relatively weak due to instrumental noise, cosmic variance and resid- 
ual foreground noise. In addition, the constraining power is further restricted by the uncertainty 
in the optical depth r to the last scattering surface. We will therefore not use the WMAP data 
in the following likelihood analysis. The BICEP data probes intermediate scales (21 < I < 335) 
around the recombination bump in the primordial .B-mode spectrum. On these scales, the pri- 
mordial signal is less affected by cosmic variance and is comparable to or larger than the lensing 
signal for tensor-to-scalar ratios r > 0.01. The QUaD experiment, whose primary aim was a high 
resolution measurement of the -E-mode signal, probes small scales (164 < I < 2026). Its ability to 
constrain the primordial signal is thus severely restricted due to lensing confusion and the rapid 
decline in the primordial signal with inverse scale. It may however be useful for constraining the 
cosmic string signal which peaks on small scales. 

The constraints plotted in Fig. Q]are all consistent with zero signal at the ~ 2a level. Detecting 
-B-mode polarization therefore remains an outstanding experimental challenge, and represents a key 
goa l for current and future CMB experiments including ground-based (BICEP-II [261 ]. QUIJOTE 



2J, PolarBear [28|, QUIET [29]), balloon-borne (EBEX [30(, Spider [3j, PIPER) and satellite 



(Planck 32|, B-Pol 33J], litebird 34|, CMBPol 35J) experiments. In what follows, we will forecast 



the constraints potentially achievable with the following five representative experiments: the Planck 
and CMBPol satellite missions, the ground-based PolarBear and QUIET (Phase II) experiments, 
and the balloon-borne experiment, Spider. In addition, for reference, we shall consider the ideal 
(but unrealistic) case where there is no foreground contamination and no instrumental noise, and 
where the lensing signal can be cleaned to around 1 part in 40 [36] . The instrumental specifications 
which we use to model the various experiments are listed in Appendix [A] Fig. [2] shows the noise 
levels of these experiments compared to the SFI signal for r = 0.03. 

C. Data Analysis Methodology 

In this subsection, we describe the methodology we use to constrain the models using current 
data, and to forecast constraints for future experiments. The parameters of the standard ACDM 
model have already been tightly constrained by CMB TT, EE and TE data [3|] and the remaining 
uncertainties in these parameters have little impact on the l?-mode power spectrum, e.g. [4]. 
Therefore, consistent with the approach adopted by Ref. [6], throughout this paper, we only vary 
those parameters which influence the level of primordial S-modes. We fix the other cosmological 
parameters at their WMAP 7-year best-fit values, which are derived under the assumption r = 
and constant n s across all wavelengths [3j: O^/i 2 = 0.02258, J7 c /i 2 = 0.1109, n s = 0.963, A s (ko) = 
2 A3 x 10~ 9 (pivot scale k = 0.002 Mpc -1 ), h = 0.71, and r = 0.088. 

1. x 2 analysis and hyper-parameters 

To constrain the three models using current data, we initially employ a conventional x 2 analysis 



to obtain the likelihood function for each data set. For LQG and SFI, we use the CAMB code 38] 
to output the transfer function for the -B-mode power spectrum. Cf can then be calculated as 

° l "4 

where Af(k) is the transfer function for each multipole I, and Pt(k) is the primordial tensor power 
spectrum (see Eqs. (fT7|) and (|3"7|) ). For the cosmic string model, we use the publicly available code 



J P t (k)Af(k) 2 dlnk, (6) 



CMBACT 



39j to generate the .B-mode power spectrum. 



We then follow the pipelines in Refs. [5|, |6] to construct the expected bandpowers for each model, 



and we use the lognormal approximation (as illustrated in e.g. [6]) to calculate x 2 according to 



xV) 



tBB 



tBB 



(a) 



\BB, 



tBB 



tBB 



(a) 



(7) 



D oa (a) 

where a is the parameter we wish to constrain, and Z BB (a) and Z BB are the model and observed 
band powers, transformed to the lognormal basis. T) BB (a) is the covariance matrix of the observed 
bandpowers, once again transformed to the lognormal basis. Minimizing the x 2 across all sampled 
values of a yields the best-fit model. 

To obtain joint constraints from BICEP and QUaD, we can simply add the x 2 values and 
minimize the resulting joint x 2 , 

Xtot = XBICEP + XQUaD- (8) 

The goodness of fit for each model can be ascertained by comparing the minimum \ 2 value with 
the number of degrees of freedom n. If the value of X 2 n in/ n is c l° se to unity within the range 



(1 — y/2/n,l + -y/2/n), we can say that the model provides a good fit to the data. If Xmin/ n -^ 
1 + \/2/n, then the model is not a good fit to the data, while if X 2 n inl n ^ 1 — y2/n, then the 
model is overfitting the data which may happen if the model has redundant free parameters and/or 
the errors on the data have been overestimated. 

Note that in the conventional joint x 2 analysis of Eq. (jSj), we have weighted each data set equally. 
This may be problematic if the two data sets are not mutually consistent, or if there are unquantified 



systematics in the data 41 



42] . In such cases, one may wish to weight the data appropriately. 



The assignment of weights often occurs when two or more of the data sets are inconsistent, and 



is usually made in a somewhat ad-hoc manner 



42] . Generally speaking, assigning the weights for 



each data set is a somewhat subjective way of performing a joint analysis, but one well-motivated 
approach to assigning weights is the "hyper-parameter" approach, formulated within a Bayesian 
context, which can objectively allow the statistical properties of each data set to determine its own 



weight in the analysis 



41 



42]. 



In the hyper-parameter technique, the effective x 2 is defined as 

^2 n i ln X 2 , 



Ahyper 



(9) 



where j sums over all of the data sets, x] is the x 2 f° r each data set, and rij is the number of 
degrees of freedom for each data set [85] . 

Once the x 2 values for the combined data set have been obtained, we can find the posterior 
distribution for the parameter a using 

-2lnP(a\D 1 ,---,D N )= X 2 , (10) 



ypcr 



4l|. 



where x 2 can be either the conventional x 2 t or the hyper-parameter version \\ 

In Appendix [Bl we calculate the expectation value and variance (see Eq. (|B17p ) of Xh yP er- This 
calculation shows that in the hyper-parameter case, a model can be said to be a good fit to the 
data if the minimum \ 2 is within the range (1 — \/V(n)/E(n), 1 + y/V(n)/E(n)), where n is the 



number of degrees of freedom for the data sets 



86]. 



2. Fisher information matrix 



In order to make forecasts for the constraints achievable with future experiments, one can 
make use of the Fisher information matrix under the assumption that each parameter is Gaussian- 
distributed. The standard Fisher matrix F a p is defined as 



43 



44] 



1. 



i;,;--Tr[C jQ C- 1 C^C- 1 ] 



"a/3 = - l-n< ,nl '■ ..;<■ I- (11) 

where C is the total covariance matrix, which includes both signal and noise contributions: 



iBB 



BB 



Cl imi l 2 rn 2 — (Cl lyS ig + N lltot )5i 1 i 2 5 mi 



7712 ' 



(12) 



Here, NPJ° t is the total noise contribution to the covariance matrix, which includes instrumental 
noise, foreground contamination as an effective noise, and confusion noise from lensing l?-modes 
(see Appendix lAl for the details). In our case where we consider only £>-mode polarization, the 
Fisher matrix can be simplified as [44j, |45(] 



F 



a/3 



£ 



2Z + 1 



sky 



(cj 



BB i j\fBB\2 



(13) 



; v / \~l,stg ' ' I, tot) 

where / s k y is the fraction of sky observed. For Planck, CMBPol, Spider and the ideal experiment, 
since these are nearly full-sky observations, we perform the summation in Eq. (|13j) from I = 2 to 
I = 3000. For the ground-based PolarBear and QUIET experiments, the summation is performed 
from I = 21 to I = 3000. We restrict the summation for these experiments to I > 20 since ground- 
based experiments are insensitive to the largest angular scales because of their finite survey areas 
(see 46J for instance). 

The inverse of the Fisher matrix F^ 1 can, crudely speaking, be considered the best achiev- 
able covariance matrix for the parameters given the experimental specification. The Cramer-Rao 
inequality means that no unbiased method can measure the i parameter with an uncertainty 
(standard deviation) less than 1/y/Fji 



43, 



441 ] . If the other parameters are not 



mown but are also 



estimated from the data, the minimum standard deviation rises to (F 1 ) 



-1^1/2 
ii 



M, 



44J . Therefore 



we can estimate the best prospective signal-to-noise ratio as a/Aa, where Aa = (F 1 )da • This 
formula will be used frequently in the following discussion. 



III. CONSTRAINING THE SFI MODEL 

A. Scalar and tensor primordial power spectra in the SFI model 

In addition to nearly scale- invariant scalar perturbations, inflationary models also predict vector 
and tensor perturbations [471 ] . However, the vector perturbations are expected to be negligible since 
these modes decayed very rapidly once they entered the Hubble horizon. We will therefore ignore 
any vector component in what follows. 

We will work in the perturbed Friedmann-Lemaitre-Robertson- Walker Universe, for which the 
metric can be written as 

ds 2 = -c 2 dt 2 + a 2 (£)(<% + h ij )dx i dx j . (14) 

The tensor perturbations h^ are described by two transverse-traceless components. The power 
spectrum for the two polarization modes of hij {hij = h^ef- + h x ef-, h = h + = h x ) is 



?TT Z 

(h k hv) = (27r) 3 <5 3 (k - k>)—P t (k), (15) 



r 2 
k 3 

where h k is the Fourier component of the perturbation field h. Standard inflationary models 
predict a nearly scale-invariant tensor power spectrum Pt(k). In order to describe the weak scale- 
dependence of Pt(k), we can define the tensor spectral index m in the usual way: 

_ dlnP t (k) 

nt = -JUT- (16) 

The tensor power spectrum can then be written in the following power-law form 

/ k \ nt 
P t (k) = A t (k ) [-) . (17) 

In the case where inflation is driven by a single inflaton field, the following calculation yields 



the primordial power spectrum of the scalar perturbations for slow-roll inflation (see 48|, |50J for 
instance. For alternative calculations, see [49].) 



P s {k) 



k=aH 



(27T0) 

9 1 V 3 



k=aH 



(2vr) 2 (3M2j)3 V' 2 

3^ \y^M&) ' 



10 



while the power spectrum of tensor perturbations is given by 



™-4(#) 



2 V 



3 7T 2 Mjj 



fc=afl" 

(19) 



Here, V(</>) is the inflaton potential, and H is the Hubble parameter at the time of inflation. It is 
customary to define the tensor-to-scalar ratio r as 

"T.-^ivf- (20) 

The tensor-to-scalar ratio and the tensor spectral index are related to the slow-roll parameter e via 



48, 



m 



r = 16c, n t = -2e. (21) 

These expressions lead to the so-called consistency relation for single field slow-roll inflation 511 ]: 



nt = ™ (22) 

Unfortunately, this consistency relation is extremely difficult to constrain observationally because 

of the small amplitude of the tensor power spectrum. We discuss the possibilities for testing this 

relation with future observations in Section IIIICI 

The normalization of the power spectrum of scalar perturbations (defined at the pivot wavenum- 

ber k = 0.002 Mpc" 1 ) is P s (k ) = (2.43 ± 0.11) x 10- 9 (la CL, WMAP 7-year data [3]). We can 

use this normalization together with Eq. (|18j) to derive the relationship between the energy scale 

of inflation and the value of r: 

i 
V\ = 1.06 x 10 16 GeV (— ) 3 . (23) 

That is, a detection of the tensor-to-scalar ratio at r « 0.01 or greater would indicate that inflation 
happened at an energy scale comparable to the Grand Unification Theory (GUT) energy scale. 

We can also use the slow-roll approximation to derive the following relation, which characterizes 
the distance in the field space from the end of inflation to the time when CMB scale fluctuations 
were created, namely the Lyth bound 



52 



53] 



M.>(JL)*. (24) 

M pl ~ V0.01/ v ; 

Thus, a tensor-to-scalar ratio greater than ~ 0.01, would directly indicate a super-Planckian field 
evolving from 4>cmb to 4> en d- Such a detection could provide important observational clues about 
the nature of quantum gravity. The boundary r ~ 0.01 is therefore an important benchmark which 
can confirm or rule out a wide class of large field inflation models. 

11 




FIG. 3: The current constraints on r from BICEP and QUaD data. 
B. Constraints on SFI from current data 

In this subsection, we present constraints on the tensor-to-scalar ratio r from BICEP and QUaD 
data. To obtain the constraints, we follow the methodology outlined in Section [II C II 

The results of the analysis are shown in Fig. [3j It is clear from the figure that the BICEP 
data provides a fairly strong upper limit on the value of r. The best-fit value is r = 0.0ll ' 2 6 (lo" 
CL), which is very close to the result obtained by the BICEP team themselves [a]. As expected, 
r is essentially unconstrained by QUaD, whose measurements are made at much smaller scales 
(I > 200) than the scale at which the primordial .B-mode signal peaks (I ~ 80). In producing 
joint constraints, we find that both the conventional x 2 an d the hyper-parameter version (Xhvper) 
are completely dominated by the BICEP data. For the conventional \ 2 analysis, the best-fit of 
the joint analysis gives r = 0.03^ ' 27 (l<x CL), and for the hyper-parameter x 2 analysis, we obtain 
r = 0.02^ ' 26 (1(7 CL). The details of the constraints are listed in the first row of Table HI Note that 
the tendency for the QUaD data to prefer larger values of r appears to be related to a marginal 
excess of power in the QUaD measurements over the multipole range 300 < I < 500 (see Fig. [1]) . 
However, the shape of this apparent "signal" is inconsistent with the SFI model (and indeed with 
any other primordial B-mode model) which suggests that it is likely due to unquantified residual 
systematics in the data rather than due to a true cosmological signal. 



Comparing the BICEP result to the WMAP 7-year results 



the tightest upper-bound on 



r that the WMAP team quote is r < 0.24 (2a). This constraint is derived from a combination 



12 



Models 


Sampling 
Range 


Conventional x 2 


hyper-parameter \ 
BICEP+QUaD 


BICEP 


QUaD (tot) 


BICEP+QUaD 


SFI: r 


(-1.0,10.0) 


n n1 +0.31+0.68 
u - ui -0.26-0.49 


10-0-3.0-9.70 


n n o+0.32+0.70 
U - UO -0.27-0.50 


n n9 +0.31+0.75 
u - uz -0.26-0.51 


LQG: m [lCT 8 M pl ] 
fc* [10- 4 Mpc- 1 ] 


(0.01, 10 2 ) 


018 + 1.16+5.74 


47.5_ 27 7 x 


022 +1.14+5.72 


20 +1.16+5.96 


(0.1, 10 2 ) 


i Q7+1- 36+5.88 


cq ,, + 131.63 x 
00 -*-23.3 X 


L07 +1.37+5.89 


L07 +1.36+5.98 


Cosmic String 
Gn x 10 7 


(10" 3 ,10 2 ) 


(o.ooi)J 5586+9 - 961 


tot- 7 fio +1 - 38+2 - 63 

lOl. 1 .VV_ 1 5 6 -3 g 






(1>500): 4.32+ 241+4 - 25 


298 +2.82+4.74 


9 on+3.45+5.69 



TABLE I: The best-fit values for the parameters, and the la and 2a CL for the scalar field inflation (SFI), 
Loop Quantum Gravity (LQG) and cosmic string models. For the SFI and LQG models, we use all the 
QUaD data and combine these with BICEP using both a conventional joint % 2 and the hyper-parameter 
approach (Xhypor)- For the cosmic strings model, in addition to the entire QUaD data set, we also examine 
the effect of removing the I < 500 QUaD data points. When combining the QUaD and BICEP data for the 
cosmic strings model, we also restrict the QUaD data to / > 500 (see the discussion in section IV Bj) . The 
notation "x" indicates that the values of the parameters are out of sampling range. 



of the WMAP data with both large scale structure measurements and the HST key project con- 
straint on the Hubble constant. It is clear that measurements of the TT and TE CMB spectra in 
combination with other astrophysical probes currently play a significant role in constraining the 
value of r. However, we note that the constraints obtained from .B-mode polarization alone are 
already comparable to the combined constraints from all other cosmological probes and are likely 
to overtake them with the next generation of CMB polarization experiments. 

In Table |H] we quote the goodness-of-fits for the various analyses and we quote the weights for 
the hyper-parameter analysis in Table Hill In Table HH E(n) is the expectation value for each fit, 
calculated using Eqs. (JB4]) and (1B17J) . If the model provides a good fit to the data, the minimum 



X 2 over the expectation value should be well within the range (l-yTfnj/Efn), 1 + y / V(n)/E(n)). 
Examining the table, we see that the SFI model can fit the BICEP data well, but that the fit to 
the QUaD data is relatively poor. This poor fit to the QUaD data adds further weight to our 
conclusion above regarding the anomalous power in the QUaD results in the range 300 < / < 500. 



C. Prospects for future observations 

In this section, we discuss the detection capabilities of future CMB experiments. In order to 
forecast the error bars of the parameters r and nt in the fiducial models, we use the Fisher matrix 
technique, introduced in Section Hi C2i 
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Goodness 
of fits 


Conventional x 2 


hyper-parameter x 2 


BfCEP 


QUaD (tot) 


BICEP+QUaD 


BICEP+QUaD 


SFI: E{n) 

X 2 mm /E(n) 
Good-fits range 


8 


22 


31 


82.58 


1.00 


1.64 


1.56 


1.26 


(0.5,1.5) 


(0.70,1.30) 


(0.75,1.25) 


(0.89,1.11) 


LQG: E(n) 

X 2 rmn/E{n) 

Good-fits range 


7 


21 


30 


75.49 


1.10 


1.76 


1.60 


1.22 


(0.47,1.54) 


(0.69,1.31) 


(0.74,1.26) 


(0.90,1.10) 


CosStr: E(n) 

xlmJEin) 
Good-fits range 


8 


18 (I >500) 


27 (QUaD I >500) 


66.6 (QUaD / >500) 


1.0 


1.35 


1.21 


1.08 


(0.5,1.5) 


(0.67,1.33) 


(0.73,1.27) 


(0.90,1.11) 



TABLE II: Reduced % 2 as an indication of the goodness-of-fit for each analysis. E(n) is the expectation 
value for each fit. If the model provides a good fit to the data, the values of xLn/^W should be within 



the range (1 - y/V(nj/E(n), 1 + y/V(n)/E(n)). 



a c ff = tia/Xa 


BICEP 


QUaD 


SFI 


1.13 


0.64 


LQG 


1.17 


0.62 


CosStr 


1.13 


0.95 



TABLE III: The value of the effective hyper-parameters a c s — tla/xa f° r the BICEP and QUaD data which 
reflect the relative weights assigned to each data set. 

In Fig. U we plot the signal-to-noise ratio (r/Ar) for a detection of tensors as a function of 
the fiducial value of r, for a number of current and forthcoming experiments. In the left panel, 
we only consider r as the free parameter, and keep nt fixed at nt = 0. We see that the Planck 
satellite can potentially detect the signal of the tensor perturbations at more than 3a confidence 
level if r > 0.05. For r = 0.1, the value of r/Ar becomes 5 which would constitute a robust 



54, 



55J. The predicted constraints 



detection. These results are consistent with those presented in 
for PolarBear, QUIET and Spider are somewhat tighter with r/Ar > 3 for models with r > 0.02. 
The predicted constraints for the proposed CMBPol mission suggest that tensor perturbations 
could be detected (at the 3<r level) for values of r as low as r ~ 0.002. Such a measurement would 
provide an excellent opportunity to differentiate between various inflationary models. We also find 
that for the ideal CMB experiment which includes only a residual lensing noise contribution (after 
de-lensing), the SFI primordial signal could be detected (at > 3<r) only if r > 10 is satisfied. 
In Fig. we also plot the signal-to-noise ratio for the combination of Planck with the ground- 
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FIG. 4: The signal-to-noise ratio r/Ar for different experiments, calculated using the Fisher matrix of Eq. 
(|13p . Left: The parameter r is treated as a free parameter but nt is kept fixed at its fiducial value; Right: 
Both r and n t are treated as free parameters. 
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FIG. 5: The signal-to-noise ratio r/Ar for different experiments, and the combination of Planck and ground- 
based experiments, calculated using the Fisher matrix of Eq. (fi"3)) . Left: The parameter r is treated as a 
free parameter but nt is kept fixed at its fiducial value; Right: Both r and nt are treated as free parameters. 



based experiments (PolarBear and QUIET). The former is sensitive to the i?-mode signal at the 
lowest multipoles £ < 20, while the latter are sensitive to the recombination peak of CP B at 
I ~ 80. Similar to 



], we find that the combination of these experiments yields little formal 
improvement in the signal-to-noise of the detection compared with the capabilities of the ground- 
based experiments on their own. However, a detection of both the recombination bump (e.g. from 
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FIG. 6: Forecasted joint constraints on the parameters of the SFI model. The contours indicate the 0070 
(la) confidence levels. The input models are indicated by the black points (Left: r = 0.01 and n t = 0, 
Right: r = 0.001 and n t = 0). 
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FIG. 7: Forecasted uncertainty on the tensor spectral index, Ant, for different experiments, calculated using 
the Fisher matrix of Eq. (fT3j) . Here, both r and nt are treated as free parameters. 



ground-based experiments) and the reionization bump (e.g. from Planck) would constitute much 
more compelling evidence for tensors than either detection would constitute on its own. 

In the right hand panels of Fig. H] and Fig. [5j we have plotted the results for the case where we 
treat both r and nj as free parameters. Comparing with the corresponding results in the left panels, 
we find (in agreement with previous works, e.g. [3a . 1571 ]) that the signal-to-noise ratios become 
much smaller due to correlations between the parameters. These correlations were investigated in 
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some detail by [56J who also explored the optimal choice of pivot scale for which the two parameters 
become decorrelated. 

In Fig. we plot the values of Ant as a function of r for various cases. Here we find that, for 
Planck, PolarBear, QUIET and Spider, the constraint on m is relatively weak unless the value of 
r is very large. For example, the predicted constraint for Planck is Ant = 0.13, and for QUIET is 
Ant = 0.18 for a model with r = 0.1. The combination of Planck and QUIET could in principle 
do somewhat better with Ant = 0.08. The proposed CMBPol mission could achieve Ant = 0.04 
while the limiting value of the ideal experiment is Ant = 0.007, which is comparable to the current 

n 

constraint on the scalar spectral index n s [3]. For lower values of r, the predicted constraints are 
correspondingly weaker. For example, for the CMBPol mission, the value of Ant increases from 
0.04 to 0.1 if we replace the r = 0.1 model with r = 0.01. 

Testing the consistency relation rif = — r/8 (see Eq. (j22j) ) is potentially one of the most powerful 
ways to test the general SFI scenario. To assess whether future experiments might achieve this 
goal, in Fig. [7] we compare the values of Ant with r/8. If Arif < r /8, then the constraint on nt 
is tight enough to allow the consistency relation to be tested. We find that Ant < r/8 is satisfied 
only if r > 0.23 for the CMBPol experiment, and only if r > 0.06 for the ideal experiment. An 
observational confirmation of the consistency relation is therefore extremely unlikely to be achieved 
with any of the currently envisaged future experiments. 



D. Single-field slow-roll inflationary models 

n this subsection, we discuss four types of single field slow-roll inflationary models, categorized 



by 



581 ]. and the implications for these models of the results presented in the previous section. 



59] 



In single field slow-roll inflationary models, the observables depend on three slow-roll parameters 



-£(£)'■ *■■*(£). *-*(^). <» 

where V(4>) is the inflationary potential, and the prime denotes derivatives with respect to the field 
<j). Here, ey quantifies the "steepness" of the slope of the potential, rjy measures the "curvature" of 
the potential, and £y quantifies the "jerk". Since the potential is fairly flat in the slow-roll inflation 
models, these three parameters must be much smaller than unity for inflation to occur. One of 
the important predictions of SFI models is that the scalar perturbations are nearly scale- invariant, 
which has already been confirmed by WMAP results 
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In SFI models, a standard slow-roll analysis yields the following relations 

r 8 16 

n t = --, n s = l + 2rjy -Gey, r = -(1 - n s ) + — rjy, a s = -24e v + 16e v rjy - 2£y, (26) 
o o o 

where n s is the tilt of primordial scalar power spectrum, and a s = dn s /dlnk is the "running" of 

n s . These formulae relate the tensor parameters nt and r to the scalar parameters n s and a s ; the 

latter can be constrained through CMB and large scale structure observations. As shown in Eq. 

(|26p . the relation between r and n s involves the slow-roll parameter rjy which in turn depends on 

the specific inflationary potential. 

The strength of the primordial tensor perturbations depends on the value of r. Observations 
have yielded quite tight constraints on n s , but we currently only have upper limits on the value 
of r. The relation between n s and r depends on the specific inflationary model, and different 
models predict very different values for r. In the following discussion, we categorize SFI models 
into four classes based on different regimes for the curvature of the potential V(4>), and discuss 
their individual constraints. 

Case A: negative curvature models rjy < 

The negative r\y models arise from a potential of spontaneous symmetry breaking. One type of 
often-discussed potentials is the form V = A 4 [1 — ((j)//j,) p ] (p > 2). This type of model predicts a 
red tilt in the scalar spectrum (n s < 1), which is consistent with the WMAP 7- year results 






addition, these models predict relatively small values for r. For the model with p = 2 in Ref. 581 ] . 

r ~ 8(1 - n s )e- N ^- ns) , (27) 

where N e is the number of e-folds, taken to be in the range A^ £ [40, 70] based on current 



observations of the CMB 3|, 160J. Here we choose the value N e = 60. Using the result n s 
0.963±0.012 [3J yields the constraint r G [0.021, 0.045]. From Fig. H (right panel), we see that this 
is close to or even beyond the sensitivity range of the Planck satellite, but is within the sensitivity 
ranges of PolarBear, QUIET, Spider and CMBPol. In other models with p > 2, the predicted 
values of r are much smaller than that of the model with p = 2. 

Case B: small positive curvature models < rjy < 2ey 

Two example potentials in this case are the monomial potentials V = A 4 (^//i) p with p > 2 for 
< rjy < 2ey and the exponential potential V = A 4 exp(^///) for rjy = 2ey. In these models, to 
first order in slow roll, the scalar index is always red n s < 1 and the following constraint on r is 
satisfied 

|(l-n s )<r<8(l-n s ). (28) 
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Using the result n s = 0.963 [3] , one finds that r £ [0.1, 0.3], which is within the sensitivity range 
of the Planck satellite, as well as that of forthcoming CMB experiments. Thus, the Planck results 
may provide some constraints on these type of models. 

Case C: intermediate positive curvature models 2ey < r\y < 2>ey 

The supergravity-motivated hybrid models have a potential of the form V — 
A 4 [l + aln(cf)/Q) + A(<^>///) 4 ] , up to one-loop correction during inflation. In this case, 

n,<l, r > 8(1 - n.) , (29) 



are satisfied. Using the result n s = 0.963 [3J, one finds that r > 0.3, which is slightly in conflict 
with the current upper limit r < 0.24 (WMAP+BAO+F , 2a CL) [3]. Fig. H shows that this 
model is also in the sensitivity range of the Planck satellite. 

Case D: large positive curvature models rjy > 3ey 

This class of models has a typical monomial potential similar to those of Case A, but with a 
plus sign for the term ((j)//j,) p : V = A 4 [1 + (</>/^) p ]. This enables inflation to occur for small values 
of <f> < M p i. The model predicts a blue tilt in the scalar power spectrum n s > 1 (Eq. (j26|) ). which 
is in conflict with current constraints on n s unless a running of n s is allowed [3] . When a running 
in n s is included, the WMAP 7-year results suggest a blue power spectrum (n s = 1.008 ±0.042 for 
1<7 CL). Therefore, even though this model is not favoured by the WMAP 7- year results for the 
case of constant n s , it is not excluded when a running of the spectral index is included. Planck 
and future CMB experiments should constrain both n s and a s to high precision and so should be 
able to definitively rule out this model. 

IV. LOOP QUANTUM GRAVITY AND ITS OBSERVATIONAL PROBES 
A. Primordial tensor perturbations in LQG models 

Loop Quantum Gravity is a promising framework for constructing a quantum theory of gravity 
in theoretical physics. Based on the reformulation of General Relativity as a kind of gauge the- 
ory obtained by 161 



621 ]. LQG is now a language and a dynamical framework which leads to a 
mathematically coherent description of the physics of quantum spacetime [63]. Constraining LQG 
theories experimentally is challenging because the quantum geometrical effect can only be tested 
at very high energy scales, beyond the reach of current accelerator experiments. In this section, 
we will calculate the possible observational signature of LQG in the CMB sky, which opens a new 
window for cosmological tests of quantum gravity. 
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There are two main quantum corrections in the Hamiltonian of LQG when dealing wit 
semi-classical approach, namely holonomy corrections and "inverse volume" corrections 
The holonomy corrections lead to a dramatic modification of the Friedmann equation as 

8ttG A p 
P. 



i the 



63, 

Q 



Si- 



H* 



-P 1 



(30) 



where p is the energy density, and p c is the critical energy density, 

p c = ^fMj 1 ~ 507.49 M 4 pl . 



(31) 



Here 7 = 0.239 is the Barbero-Immirzi parameter, which is derived from the computation of the 



black hole entropy 66j. Note that we use the reduced Planck mass in our calculation. 

A generic picture for this model with the holonomy correction is the bouncing behavior exhibited 
when the energy density of the Universe approaches p c . The negative sign in Eq. (|30p is an 
appealing feature in the framework of LQG such that the repulsive quantum geometry effect 



becomes dominant in the Planck region 



63 



641 ] . This triggers a contraction period before the 



bounce, during which time the Hubble parameter is negative and the Hubble radius is shrinking. 
As a result, the perturbation modes on the largest scales crossed the Hubble horizon and froze 
out during the contracting period, until the end of the contracting stage when the Hubble horizon 
increased again. For the very large scale modes, this pre-inflationary bounce may imprint distinctive 



64]. 



features in the CMB sky, since they stretched out of the horizon at very early times 63, 

Unfortunately, the power spectrum for the scalar perturbations is somewhat hard to obtain 
because in the case of the holonomy correction, the anomaly free equations are still to be found 
671 ] . Therefore, in the following discussion, we will focus on the tensor power spectrum of LQG 
and pursue the constraints obtainable on the model from I?-mode observations. 

Due to the pre-inflationary contracting period, the tensor power spectrum for LQG can be 



calculated numerically. In 64j], a simple parameterized form of the power spectrum is introduced 
as follows 

2 



Pt 



7T- 



H 



1 



1 + 



AR-2 



i + 0/fc*) 



(32) 



1 + (h/k) 2 

where H is the Hubble parameter during the inflationary stage, £;* is the position of the highest 
peak in the power spectrum, and the quantity R is related to the mass of the scalar field as 

0.64 



R 



-in 



,0.32 



Mpl 
in 



(33) 
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It is interesting to note that Eq. (j32|) reduces to the SFI result of Eq. (fT9[) for fe* — > 0. In this 
paper, for simplicity, we consider the tensor power spectrum (Eq. I3"2"j) with a constant H in the 
early stage of inflation, which corresponds to the specific case of de Sitter inflation. 

In this model, we assume that inflation is driven by the potential V((p) = m 2 (p 2 /2. The Hubble 
parameter is related to this potential via 

1 1 



H' 



3MJ, 



v{<p) 



6M2, 



2 i.2 

m 0j, 



(34) 



where (pi is the initial value of the scalar field at the beginning of inflation. The number of e-folds 
can be calculated as (using the slow-roll approximation 2>H(p = —V) 

"f 



N r . 



Hdt 



"4 f ^ 



V 



1 



4MJ, 



?)■ 



(35) 



Since (pi S> </>/, the above equation is approximately given by cp 2 ~ 4iV e Mpj. The Hubble parameter 
is then (from Eq. (I34p ) 



H 2 ~ -N e m 2 . 
3 



(36) 



Substituting this into Eq. ([32]) . we arrive at the following expression for the tensor power spectrum: 



Pt 



AN P 



m 



1 



1 + 



4 x (Svr) -^ x (m/M p i) 



-0.64 



(37) 
to project the 



3^ V m pi7 l + (h/k) 2 
where N e is the number of e-folds which we fix at -/V e ~ 60. Finally, we use Eq. 
perturbation modes onto the CMB sphere to find Cf for the LQG model. 

In the left panel of Fig. [51 for a number of representative sets of parameters, we plot the 
primordial tensor power spectrum for the LQG model alongside the signal expected in an SFI 
model for a number of different values of r. For the SFI model, since the power spectrum tilt 
ni is very small {rit = for de Sitter inflation), the power spectrum is very flat on all scales. 
In comparison, the tensor power spectrum of LQG exhibits a bump feature, which is in fact the 
signature of the pre-inflationary contraction period. Very large scale modes stretched out of the 
Hubble horizon during the contracting period before inflation, and can be described by the solution 



in the Minkowski vacuum fk = e 
takes the form Pt{k) ~ k 3 \ fk \ 2 ~ k 2 



63 



63, 



64|. Thus, the power spectrum at very large scale 



64j. In contrast, the small scale modes are well within the 



Hubble horizon and so the power on small scales is similar to the scale-invariant power spectrum of 
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FIG. 8: Left: the primordial tensor power spectrum for different models of inflation. Right: the correspond- 
ing i?-mode angular power spectrum generated by different inflationary models. Note the lensing i?-mode 
signal which acts as an effective noise for detecting the primordial Cf B signal. The units for m and fc* are 
M p i and Mpc -1 respectively. 

the SFI model. The bump in the power spectrum on larger scales is characterized by the magnitude 
of &*. As k* increases, the bump is shifted to smaller scales and vice versa. The amplitude of the 
spectrum and the width of the bump are determined by the mass parameter m. We will link these 
two important parameters to the energy scale of inflation and the current Hubble horizon scale in 
the next subsection. 

The right panel of Fig. [8] shows that the bump in the primordial power spectrum results in a 
peak in the CMB l?-mode power spectrum, which is slightly different to that of the SFI model. In 
addition, if the peak of the LQG spectrum is normalized to the same magnitude as that of SFI, 
the small scale power will be suppressed in the LQG model, as compared to the SFI model. 

B. Constraints from current data 



In this subsection, we use the BICEP and QUaD data to constrain the parameters of LQG 
models. Before we perform the parameter estimation, we link the two parameters m and k* with 
the energy scale of inflation, and with the current Hubble horizon scale. 

The parameter m relates to the energy scale of inflation in LQG as follows 

1/2 



V~i = 3.02 x 10 15 GeV 



m 



(38) 



.lCJ- 7 M pl , 

where 10 15 GeV is around the GUT energy scale. Therefore, a detection of m > 10~ 7 M p i would 
strongly suggest that the energy scale of inflation is above the GUT scale. 
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FIG. 9: Probability distribution function (PDF) for the parameters in the LQG model. Left: constraints on 
the mass parameter rh — m/M p i. Right: constraints on k*. The red, green, and brown curves overlap with 
each other, indicating that the vast majority of the constraining power for the combined data sets comes 
from the BICEP data. 

The parameter k* describes the position of the peak in the primordial power spectrum. We 
can compare it with the current Hubble wavenumber, which is kjj = Hq ~ 2.33 x 10 Mpc -1 . If 
fc* > ksj then modes with physical wavelengths (A*) equal to the Hubble horizon at the beginning 
of inflation will have wavelengths less than the current Hubble horizon, whereas if k* < kg their 
wavelengths will be larger than the current horizon scale. Thus, if fc* > kij, we would expect to be 
able to find pre-inflationary fluctuations within our current Hubble horizon 
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641 ] . Conversely, as 



k* —7- 0, the primordial tensor power spectrum (Eqs. (|32p and (|37p) reduces to the scale-invariant 
tensor power spectrum as noted above. Therefore, a non-zero detection of k* would strongly 
indicate the existence of a bounce and of a contracting period before inflation. 

The current constraints on the parameters k* and m are shown in Fig. [9j In the left panel, we 
have marginalized over the parameter k* and we plot the PDF for the mass parameter m. The la 
upper bound is m < 1.36 x 10~ 8 M p i. The detailed results are listed in the second and third rows 
of Table [H We note that the GUT scale mass m ~ 10 -7 M p i is excluded at the la level, but is 
still well within 3a. As was the case with the SFI models, the small-scale QUaD data is unable to 
constrain the LQG models. Once again, the combined constraints are dominated by the BICEP 
data as is clear from the figure and from the results listed in Table 1. 

In the right panel of Fig. [U we show the PDF for k* (marginalized over the m parameter). 
Once again the results are dominated by BICEP. From this plot, wee see that there is a peak 
in the PDF at k* = 1.07 x 10 Mpc . Although it is not statistically significant, a detection 
of such a feature would be an interesting result for the pre-inflationary bouncing behavior, since 
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FIG. 10: 2D constraints on the parameters m = m/M p i and fc*|Mpc x ] of the LQG model. The red line 
(BICEP ler CL) overlaps with the brown line (BICEP+QUaD, la CL). 

the bounce of the primordial tensor power spectrum is characterized by a non-zero k* as we have 
already discussed. We further note that the peak value of fc* is only slightly smaller than the 
current Hubble wavenumber fc#, which would indicate that modes which had the same length 
scale as the Hubble horizon at the beginning of inflation have not evolved into the Hubble horizon 
yet. In that case, the bump in the tensor power spectrum of LQG is a super-horizon feature. 
However, we stress again that all of our results are upper limits only and our formal constraint is 
k* < 2.43 x lO^Mpc" 1 (la CL). 

In Fig. [TUl we plot the two-dimensional constraints on the parameters m and k* on a log scale. 
Clearly, the current data is unable to provide strong constraints on the joint distribution of these 
two parameters and can only provide upper limits. Once again, as expected, the constraints are 
dominated by the BICEP data. 

C. Prospects for future experiments 

We now investigate the prospects for detection of LQG signatures with future CMB experiments 
by studying the projected constraints on the two parameters m and k*. Once again, we use the 
Fisher matrix approach described in Section III C 21 and as in Section IIIICl we fix the background 
parameters at their WMAP 7- year best-fit values [3J. 

In Fig. [TT] we plot the projected signal-to-noise ratio for the parameters m (left panel) and 
fc* (right panel) for forthcoming experiments. In the left panel, we plot the signal-to- noise ratio 
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FIG. 11: Predicted signal-to-noise ratios for the parameters in LQG for forthcoming and future experiments. 
Here rh is the mass parameter, and k* is the position of the bump in the BB power spectrum. 
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FIG. 12: Forecasted 2D constraints for the LQG model. The input models are indicated with the black 
points. (Left: m = 10~ 8 and fc* = 0.002MPC- 1 . Right: m = 10~ 8 and K = 0.0002MPC- 1 ). 

m/Am as a function of m. For this plot, we have kept /c* fixed at k* = 0.002 Mpc^ 1 . We find 
that, if m > 10~ 7 M p i, i.e. the energy scale of inflation is higher than the GUT scale, then 
Planck, Spider and CMBPol could potentially detect the LQG signal at more than 10o" due to 
their large sky coverage. In contrast, the smaller scale experiments (PolarBear and QUIET) would 
only detect the signal at the 2 — 3a level. For m = 10~ 8 M p i, close to the upper bound of the 
current la confidence level, the large scale survey experiments (Planck, Spider and CMBPol) can 
still detect the signal at more than 5a. For this mass parameter and value of &;*, the ground-based 
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experiments, PolarBear and QUIET, would be insensitive to the signal since the S-mode power 
spectrum from LQG falls off extremely rapidly with increasing 1. However, for larger values of fc*, 
the peak of the LQG Cf B power spectrum moves to smaller scales. We therefore expect a general 
trend whereby the large-scale experiments will be sensitive to models with small values of k* and 
small-scale experiments will be sensitive to models with larger values of fe*. 

This is illustrated in the right-hand panel of Fig. [TT1 where we plot forecasts for A^/A/c* as a 
function of k*. For these results, we have fixed m = 10~ 8 M p i. We find that the signal-to-noise 
ratio of k* does not monotonically increase with increasing k* for any single experiment. Since 
this parameter controls the angular scale at which the LQG .B-mode signal peaks, as we vary 
k* we move between the sensitivity ranges of different experiments. For example at k* ~ 0.002 
Mpc -1 , both Planck and Spider could detect the signal at more than 5a whereas PolarBear and 
QUIET would achieve only a marginal detection. Conversely, if k* ~ 0.02 Mpc -1 , the reverse is 
true: PolarBear and QUIET would make strong (> 5a) detections while Planck and Spider would 
struggle to detect a signal. 

In Fig. [T2| we plot the two-dimensional constraints on the parameters m and k* for two typical 
models. The left panel shows the forecasted constraints for a model with k* = 0.002 Mpc -1 and the 
right panel shows the constraints for a model with k^ = 0.0002 Mpc -1 . In both cases, a fiducial 
value of m = 10~ 8 M p i was adopted. We find that the former case can be well constrained by 
either Spider, CMBPol, Planck+PolarBear or Planck+QUIET while the latter case can only be 
meaningfully constrained by the CMBPol mission. 

V. COSMIC STRINGS AND THEIR DETECTION 

A. B-mode polarization from cosmic strings 

Cosmic strings have been proposed as a possible source of the inhomogeneities in the Universe 



681 ] . Although current observations of the CMB temperature and polarization power spectra 
suggest that it is inflation, rather than cosmic strings, which is the main source of the primordial 
density perturbations [3[ , there is still significant motivation to search for the signature of cosmic 
strings from both theoretical and observational considerations. 

Cosmic strings can be formed in several inflationary pictures, and particularly in Brane inflation 



models 



iiiU, 



69( | . Brane inflation arises from the framework of high dimensional string theory 



and is a further important model for sourcing the dynamics of inflation. In this model, the high 



2G 



dimensional Brane and anti-Brane collided and annihilated and cosmic strings were produced at 
the end of the inflationary epoch. If this scenario is correct, the resulting strings would have made 



ZQfl. 



an observable imprint on the CMB sky by way of the Kaiser-Stebbins effect 

Although current observations suggest that inflation sources the majority of the CMB 
anisotropy, one cannot rule out a significant (up to ~ 10%) contribution from cosmic strings 71] ■ 



In this section, we will use CMB data to constrain the level of cosmic strings. However, in contrast 
to other works (see e.g. [Til] ) , we focus on the possible detection of cosmic strings through the 
.B-mode power spectrum alone. U-mode polarization can only be generated in the early Universe 
by vector and tensor perturbations, which provides a complementary route for detecting cosmic 
strings. Although the contribution of scalar perturbations from cosmic strings is subdominant 
(< 10%) compared to the contribution from SFI, the contributions of vector and tensor pertur- 
bations from cosmic strings may constitute a very significant fraction of the l?-mode polarization 



72]. 



power on small angular scale (high multipoles) 

To predict the Cf power spectrum generated by cosmic strings, one must understand the 
evolution of a cosmic string network and it is important to know the characteristics of the scaling 
regime which needs to be assumed in the numerical simulation. There are two popular ways of 
making progress. One is to solve the Nambu equations of motion for a string in an expanding 
universe and ignore the effects of radiation backreaction (hereafter the Nambu-String model); the 
other is to solve the equations of motion for the Abelian-Higgs (AH) model, but to limit the 
dynamical range of the simulation. In this work, we will consider only the Nambu-String model 
but we note that the Abelian-Higgs model can be constrained in a similar way. 

In order to calculate the 5-mode power spectrum, including the contributions of both vector and 



tensor perturbations, we use the publicly available code CMBACT [39J, 



401 ] to generate a fiducial 



CP B for cosmic strings with the tension of the strings set to G^lq = 10 -7 . Since the amplitude 
of the S-mode power spectrum generated by cosmic strings is simply proportional to the square 
of the cosmic string tension, we can scale the fiducial spectrum to any other value for the cosmic 
string tension using 

cr=cr fi (§£)\ (39) 



dd n 

where C l ' is the power spectrum normalized at G/iq = 10" 
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FIG. 13: Current constraints from B-mode polarization on the cosmic string tension. To obtain the con- 
straints, we have fixed the wiggling parameter at a = 1.9. 

B. Constraints from current data 

In Fig. [TBI we show the current constraints on the cosmic string tension Gfi from the BICEP and 
QUaD Cf B data. The upper bound from BICEP alone is Gfi < 9.961 x 1(T 7 (2a CL). Using the 
QUaD data alone the result is G\i = 7.60^ 3 ' 60 x 10~ 7 (2<r CL), shown as the blue curve in Fig. 1131 
Taken at face value, the QUaD result represents a 2.8a detection of the cosmic string tension. 
Referring back to Fig. [1] and the discussion in Section IIIIB| this result is likely coming from the 
apparent excess of power seen in the QUaD data on scales 300 < I < 500. To investigate further, 
we have repeated the analysis with the I < 500 QUaD data removed. This results in the likelihood 
function for Gfj, shown as the brown curve in Fig. [T3j Excluding the / < 500 QUaD data shifts the 
peak of the likelihood significantly towards a smaller value (best-fit G/j, = 4.32 x 10 -7 ), and the 
constraint is now consistent with zero at the la CL. The 2a upper limit becomes Gfi < 8.57 x 10 
(see Table [J for the full set of results). The fact that the QUaD result changes drastically when we 
remove the I < 500 measurements suggests a problem with the I < 500 data. As in the case of the 
SFI constraints presented earlier, we note that the shape of the QUaD data at I < 500 is clearly 
inconsistent with the expected .B-mode signal for cosmic strings. Once again, we suspect that the 
anomalous signal seen in the QUaD data between I ~ 300 and / ~ 500 is likely due to unquantified 
systematics. We therefore consider the QUaD result restricted to / > 500 to be a more robust 
constraint and consequently we quote this as our main result. 

We now examine the constraints obtained from combining the BICEP data with the I > 500 
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QUaD data. Since the peaks of the two likelihood functions do not overlap, we will carefully 
consider both the conventional x 2 analysis and the hyper-parameter \ 2 - When we use the conven- 
tional x 2 ) the peak of the combined likelihood lies midway between the peaks of the two individual 
likelihoods (green line in Fig. [T3j) . and the resulting constraint is G[i < 7.72 x 10 -7 (2a CL). If 
we use the hyper-parameter \ 2 •, the peak of the joint distribution moves slightly further towards 
zero, and the best-fit is Gfi < 8.01 x 10~ 7 (2a CL). Since the conventional and hyper-parameter 
approaches give very similar results, this suggests that the BICEP and I > 500 QUaD data are 
mutually consistent which adds confidence to the joint constraint. 

A number of previous works have also attempted to constrain cosmic strings throu gh their 



7MIM- In 



imprint on CMB (TT,TE,EE), large scale structure and gravitational waves data 

a recent work [71[, constraints on the cosmic string tension were obtained from a combination 

of CMB data (including the WMAP 5-year, ACBAR, BOOMERGANG, CBI, QUAD and BIMA 

observations), matter power spectrum data from the SDSS Luminous Red Galaxies sample, and 

Big Bang Nucleosynthesis constraints on the baryon fraction from measurements of deuterium at 

high redshift. They obtained a combined upper limit of Gfi < 2.2 x 10 -7 (2a) in the Nambu-String 

case. 

As expected, this limit is much tighter than what we have obtained from S-modes alone since 
current measurements of the TT, TE and EE power spectra are much stronger than the BB data 
that we have used here. However, we note that in the analysis of Ref. 71], the main constraining 



power comes from the CMB TT and SDSS data and ultimately constraints from such data will 
be limited by degeneracies with other parameters (most notably, the spectral index n s ). As is the 
case for SFI models, the advantage of using the BB power spectrum to constrain the cosmic string 
tension is that Cf B is only very weakly dependent on the other cosmological parameters, e.g. 
n s . Our constraint is therefore an independent check of other constraints obtained on Gfi using 
different data and our result, although weaker, is consistent with previous analyses. Cosmic string 
constraints from forthcoming CMB polarization experiments will likely close the gap with other 
techniques in terms of constraining power. We now turn to examining the constraints achievable 
with these forthcoming experiments. 

C. Prospects for future experiments 

In this subsection, we forecast the detectability of the cosmic string tension Gfi for future 
experiments. Once again, we use the Fisher matrix formalism described in Section III C 21 In 
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FIG. 14: Forecasts of the signal-to-noise ratio for the cosmic string tension Gfj, potentially achievable with 
future CMB polarization experiments. 

Fig. PT41 we plot the signal-to- noise ratio for a measurement of G[i for the various experiments. As 
one might expect, the Planck and Spider experiments are unable to tightly constrain the cosmic 
string tension since neither experiment will produce sensitive i3-mode measurements on small scales 
where the string signal peaks. 

For example, the Planck satellite can detect a cosmic string signal at the 3a level only if Gjjl > 
2x 10 -7 is satisfied. The ground-based experiments, PolarBear and QUIET, will obviously be better 
at constraining Gfi than Planck and Spider. PolarBear should be able to detect G\x > 5.0 x 10 -8 
(at the 3d level) while QUIET should be able to detect G/i > 4.0 x 10 -8 (again at 3cr). Adding 
Planck to either of these experiments does not change the results significantly. CMBPol is much 
more sensitive than the other experiments for all fiducial G[i values, and its 3<r detection limit is 
G\x > 1.5 x 10~ 8 . Finally, we find that the ideal CMB experiment can detect cosmic strings at the 
3<t level if G[i > 1.5 x 10~ 9 is satisfied, which represents the fundamental detection limit for CMB 
-B-mode polarization experiments. 

VI. CONCLUSION 

In this paper, we have examined the observational signatures of three different models of the 
early Universe related to the inflationary process. These three models are in turn motivated by three 
different aspects of microscopic physics: single field slow-roll inflation from effective field theory 
(SFI), loop quantum cosmology from loop quantum gravity (LQG) and cosmic strings from Brane 
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inflation and/or high dimensional string theory. We have discussed their potential observational 
signatures in the I?-mode polarization of the CMB, and we have constrained the parameters of each 
model using the latest CMB polarization data from the BICEP and QUaD experiments. Using a 
Fisher matrix formalism we have forecasted the constraints achievable on these models using future 
CMB polarization observations from a number of experiments including Planck (space), Polar Bear 
(ground), QUIET (ground), Spider (balloon), CMBPol (space) and an idealized experiment. 

We first discussed the SFI model. From the Lyth bound relation, we know that r ~ 0.01 is 
an important bound for inflation reaching the GUT scale, and the inflaton evolves over the trans- 
Planckian region. The constraints we obtained from current S-mode measurements are shown in 
Fig. El Using the BICEP data alone, we find r = O.OIIq'^ (la CL) in close agreement with the 
BICEP team's own analysis [6J. As expected, this constraint does not change significantly on the 
addition of the small-scale QUaD data. Looking to the future, we find that the Planck satellite may 
be able to detect r ~ 0.05 (at the 3a CL), while the Spider, QUIET and PolarBear experiments 
all have the potential to make a 5<r detection of r ~ 0.05. The possible future satellite mission, 
CMBPol could detect r ~ 0.01 at the 20a CL, and could even detect r ~ 0.002 at the ~ 3a CL. 
(All of these forecasts are for the case where the tensor spectral index is held fixed at n% = 0.) In 
summary, we find that all future experiments can potentially constrain the value of r with sufficient 
sensitivity to allow the SFI model from effective field theory to be tested in an interesting way. 

We have also discussed the LQG model, which predicts a pre-inflationary bouncing era. Before 
the bounce, the Universe was contracting and dominated by vacuum energy. Its tensor power 
spectrum is characterized by two parameters: m and &*, where the mass parameter m controls 
the magnitude of the B-mode power spectrum, and k* controls the scale of the peak in the B- 
mode spectrum. Our joint likelihood analysis using current data yields m < 6.16 x 10 -8 M p i and 
k* < 7.05 x 10 Mpc -1 (both 2a CL upper limits). The PDF for the parameter k* exhibits 
a peak at position k* = 1.07 x 10~ 4 Mpc -1 . Although this peak is not statistically significant, 
we note that were the value of k* to be around this value, this would constitute evidence for a 
pre-inflationary bounce. We have also presented forecasts for constraining the LQG model using 
future CMB experiments. We find that if the value of h* is as large as 0.002 Mpc -1 , a number of 
ongoing and future experiments (Planck, Spider and CMBPol) could potentially detect the signal 
as long as m > 2.5 x 10 -9 M p i. However, if the value of k* were to be as low as 0.0002 Mpc -1 (as 
is mildly indicated by current data), then the signature of LQG becomes quite difficult to detect. 
We find that, for a typical choice of m = 10 -8 M p i, only CMBPol and the ideal experiment could 
detect the signal for such a low value of A;*. 
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Finally, we have presented current and prospective constraints on the cosmic string tension. We 
find that the BICEP and QUaD data constrain the cosmic string tension to be G/j, < 9.96 x 10 -7 and 
G/i < 8.57 x 10~ 7 respectively (both 2a CL upper limits). The combined constraint is Gfj, < 8.01 x 
10~ 7 , which is weaker, but comparable to the constraints from the CMB temperature anisotropy 
power spectrum. In terms of forecasts for the future, we find that the high-resolution ground- 
based experiments, PolarBear and QUIET, are more useful for constraining Gfi (as compared to 
e.g. Planck and Spider) since they are much more sensitive to the B-mode power spectrum on 
small angular scales. These two experiments could detect G\i ~ 1CP 7 at more than 10a CL. We 
also find that CMBPol can detect the signal of cosmic strings if G/i > 1.5 x 10~ 8 at 3er CL, and 
the ideal CMB experiment can detect the signal at the 3<r level even if the tension of cosmic string 
is as low as Gfi = 1.5 x 10 -9 . 

Although the i?-mode polarization derived constraints which we have presented in this paper 
are currently only upper limits, they are nevertheless already comparable to the equivalent con- 
straints obtained using a combination of all other cosmological data. In terms of constraining the 
parameters of early Universe models, 5-mode polarization is clearly a very powerful tool and will 
likely overtake other early Universe probes with the advent of the next generation of CMB polar- 
ization experiments. In this paper, we have not directly considered the issue of model selection. 
However, it is likely that, in addition to constraining model parameters, future sensitive 5-mode 
observations will also allow us to distinguish between models of the early Universe such as those 
considered in this paper. 
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Appendix A: Instrumental Characteristics of CMB Experiments 

To calculate the total noise power spectrum N^ B (Section III C 2p . we require the experimen- 
tal specifications of each experiment, including the levels of both residual foreground noise and 
instrumental noise. We list the instrumental noise for each frequency channel the experiments we 
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Band center [GHz] 


30 


44 


70 


100 


143 


217 


353 


FWHM [arcmin] 


33 


24 


14 


10.0 


7.1 


5.0 


5.0 


N£* t (i) [10-VK 2 ] 


2683 


2753 


2764 


504 


279 


754 


6975 


/sky 


0.65 



TABLE IV: Instrumental parameters for the Planck satellite (space-based experiment) [32j. Here we have 
assumed 4 sky surveys (28 months). 



Band center [GHz] 


90 


150 


220 


FWHM [arcmin] 


6.7 


4.0 


2.7 


#&» [10-VK 2 ] 


5.2 


4.3 


44.0 


/sky 


0.024 



TABLE V: Instrumental parameters for the ground-based PolarBear experiment 



Band center [GHz] 


40 


90 


FWHM [arcmin] 


23 


10 


NfJ/i) [10-VK 2 ] 


0.26 


0.64 


/sky 


0.04 



TABLE VI: Instrumental parameters for the ground-based QUIET experiment 
the phase-2 experiment. 



29] . Here we have assumed 



Band center [GHz] 


100 


145 


225 


275 


FWHM [arcmin] 


58 


40 


26 


21 


A&?,(i) [10-VK 2 ] 


84.4 


47.4 


395 


1170 


/sky 


0.5 



TABLE VII: Instrumental parameters for the balloon-borne Spider experiment [31| . Here we have assumed 
a 30-day LDB flight. 



Band center [GHz] 


30 


45 


70 


100 


150 


220 


340 


FWHM [arcmin] 


26 


17 


11 


8 


5 


3.5 


2.3 


Nffi® [10-VK 2 ] 


31.21 


5.79 


1.48 


0.89 


0.83 


1.95 


39.46 


/sky 


0.8 



TABLE VIII: Instrumental parameters for the mid-cost (EPIC-2m) CMBPol satellite mission 35]. 
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Parameter 


Synchrotron 


Dust 


A S ,D 


4.7 x 10~ 5 /iK 2 


1.2 x 10~ 4 /jK 2 


vq 


30 GHz 


94 GHz 


h 


350 


900 


a 


-3 


2.2 


pBB 


-2.6 


-1.4 



TABLE IX: Assumptions for foregrounds parameters 



have considered in Table HVTlVIIIi The effective noise power spectrum Nj^ B is given by the optimal 



combination of the channels 

rBBl-2 



35] 



[Nf 






T BB , 



T BB 



BB , 



(wj(i) + N^mwjij) + n««(j))-(i + <%) 



(Al) 



where iVj^(i) and NPf(i) are the instrumental and residual foreground noise power spectra, 
respectively. Note that the noise power spectra iVjj^(i) listed in the tables do not include the 
window function of the instrumental beam exp[l(l + 1)6^/(8 ln2)]. 

To model polarized foregrounds, we focus on diffuse synchrotron (S) and dust (D) emission. 
The foreground contamination can be quantified by the parameter a g which multiplies the power 
spectra Cgf(i), C^(i) of the foreground models. The smaller the value of <r fg the deeper the 
foreground cleaning. Throughout this paper, we adopt a s = 0.1. The residual foreground noise is 
given by (see [35] for instance) 



*£?(0 



£ [cff®o*+tf?f(ij 



f=S,D 



(A2) 



where Nfi(i) is the noise power spectrum arising from the cleaning procedure itself in the presence 
of instrumental noise. 
Following 
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79, 



and dust emission as 



80J, we model the scale (/) and frequency (vi) dependence of the synchrotron 

2as / ; \ As 



cg?(0 



A s 



Vi 



V{) 



Pi 



(A3) 



and 



cEH(i)= P 2 A D 



V; 



"0 



2a L 



l\ P ° 



e hu /kT _ 1 
phui/kT ^ 



(A4) 
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In Eq. (|A4j) . p is the dust polarization fraction, estimated to be 5% [79(, and T is the temperature 
of the dust grains, assumed to be constant across the sky with T = 18K [79(]. Other parameters in 



Eqs. (|A3p . (|A4p are specified in Table HXl taken from 35] . 

The noise term N?P{i) (/ = S, D) entering Eq. (|A2|) is calculated in [3a. 17' 

ra chan (n ch an - l)/4 V^ref / 

Here, n c h a n is the total number of frequency channels used in making the foreground template 
map, and z/ re f is the frequency of the reference channel. In the case of dust, i/ re f is the highest 
frequency channel included in the template making, while in the case of synchrotron, z/ re f is the 
lowest frequency channel. The value of a is given in Table IIXI for different foreground models. 
We note that the ground-based experiments are insensitive to the largest angular scales, so when 
calculating the Fisher matrix using Eq. (|13|) . we sum over the I from 21 to 3000. In addition, 
for ground-based experiments, the small scale fluctuations are not very sensitive to the residual 
foreground noise, and we can also pick out relatively clean patchs of sky where the foreground 
contamination is minimal. Therefore, to forecast the results for PolarBear and QUIET, we have 
not included a residual foreground noise term. 

In addition to instrumental and residual foreground noise, gravitational lensing converts E- 



mode go 
signal 



arization into l?-modes on small angular scales, contaminating the primordial i3-mode 



811 ] . The lensed C, (lens) will also contribute to the total noise power spectrum TV, . The 
total noise power spectrum therefore becomes 

Nfi* = N? B + Cl iE '(lens). (A6) 

For the ideal case, we assume that there is no instrumental or foreground noise, and that we 



can successfully de-lens the CMB observations to a level of about 1/40 of the lensing signal 361 ] . 
In this case, the total effective noise power spectrum is 

Nf™ = 1/40 x Cf B (lens). (A7) 

Finally, we adopt / s k y = 0.8 for the ideal experiment, which is the same as that used to model 
CMBPol. 

Appendix B: Statistics of the conventional \ 2 and the hyper-parameter x 2 

In this appendix, we first review the basic results of conventional x 2 statistics and then generalize 
the analysis to the hyper-parameter technique. 
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1. x 2 statistics 

A conventional joint x 2 analysis will minimize the following combined x 2 



xL = X>i' (Bl) 

3 

where each Xn follows the chi-square distribution 

f(xl) = -K^ixD^eM—xD- (B2) 

It is easy to show that this chi-square distribution is properly normalized, i.e. 



f(xl)dxl = 1, (B3) 



/o 
and the expectation value and the variance are 



E(x 2 n ) = n, V{x 2 n ) = In. (B4) 

Therefore, the minimum x 2 value of a properly constrained model should satisfy the following 
relation 



1 _v^R < ^L^ <1+ v^R , B5) 

E(n) ~ E(xl) ~ + E(n) " (B5) 



For the x 2 with order n, this is 

1_ Jl < ^H«n <l + J I, (B6) 

V n n V n 

If the Xmin/ n — 1 + J^,we can say that the model does not provide a good fit to the data, whereas 
if Xmin/ n — 1 — \/n' we sa y that the model overfits the data, which may mean that the model has 
redundant free parameters. 

If there are m constraints on the n random variables, then Xn still follows the chi-square distri- 
bution, but with order n — m [82] 

fixl) = n- m l — ^ (Xn)^' 1 *M-\xl)- (B7) 

2^r(M) 2 

It is straightforward to verify, that the shape of the distribution does not change, but the expec- 
tation value and the variance are changed simply as n — > n — m. 
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2. Hyper-parameter x 2 

The hyper-parameter approach to combining the constraints from different data sets can be 
useful in the case where the different data sets have different levels of systematics. To weight each 
data set, one should multiply each x 2 by a f ree parameter, 

xlyper = J2 a jX 2 j, (B8) 

3 

where aj is the weight parameter for each data set. One can marginalize these weight parameters 
in a Bayesian analysis, and in [41|, the authors found that instead of minimizing the combined x 2 > 
one should instead minimize the following quantity 

Xhyper = ^ n J ln 4' ( B9 ) 

3 

where n,- is the number of degrees of freedom for each data set. If we only consider one data set, 
the hyper-parameter statistic becomes 

Xhyper = nm X^- (BIO) 

From Eq. (|B2|h using the following transformation 

fy = fx 



dY 



dX 
one can show that the distribution of the hyper-parameter statistic is given by 



(Bll) 



dixlypev) = n ,„, exp(-Xhyper) ex P(- ~ exp(-Xh yp er))> ( B12 ) 



which has already been properly normalized. (When calculating the integral, one should use the 
transformation exp(^) = y). One can then show that the expectation value and variance of the 
Hyper-parameter distribution is 



E(xL ve r)=n(ln2 + i; (^)), V( X L per ) = n^i ( ~) , (B13) 



where ip n (x) is the "digamma function" defined as derivatives of the log Gamma function 

jn+l 
^ = d^ lUr{X) - (BW) 

If there are m constraints on the n random variables (e.g. m parameters), then the distribution 
of the conventional Xn follows Eq. (|B7j) . It is then easy to show that the form of the distribution 
is unchanged only if the Xh yP er ^ s defined as 

Xhyper = ( n ~ m ) ln xl- (B15) 
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Proof From Eqs. (|B7|) and (jBlip . one can derive the following PDF 

f(xl) 



9 VAhyper; 



'd X 2 



hyper 
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\ V n — m I I n — m \n — m 

1 
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exp 



ypcr 



n — m 



exp 



A. hyper 

n — m 

1 



n — ml o^^r ( ra ~ m ) 
i 



n — m _ 
2 



f 1 / ^hy 

exp --exp 

\ 2 \ n — m 



Hyper 



J- ■!■ / A.hyper \ / J- / Xhyper 

n-m 2 ^ rfM) eXP I -2~ J 6XP "2 eXP U 3 ^ 



(B16) 



Therefore, in the case of m constraints, the distribution keeps its form, and the expectation value 
and variance become 



^(Xhyper) = (n - m) [ lli 2 + 1p 



. ^(Xhyper) = (n - m) Vl 



(B17) 



Thus, to ensure that the form of the Xh yP er distribution function is unchanged, the xL per needs to 



be defined as 



Ahyper 



ypcr 



^njhixlp 



(B18) 



where rij = nd ata — m is the number of degree of freedom. 

The hyper- parameter approach is an objective way to weight each data set when producing 
joint constraints. The value of the weight is simply the value of the effective hyper-parameter, 
which is defined as 
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a a = — 9", 

xi 



(B19) 



where A specifies a particular data set. Therefore, the larger the value of a, the larger the weight 
that the particular data set takes (see Table IIIip . 
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